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ABSTRACT 
Forces  Due  to  Waves  on  Submerged  Structures.   (May  1970) 
George  Edward  Shank,  B.S.C.E.,  Duke  University 
Directed  by:  Dr.  John  B.  Herbich 

Forces  due  to  gravity  waves  on  several  models  of  half- 
cylindrical  and  rectangular  shape  were  determined  experimentally 
in  a  two  dimensional  wave  tank.  The  parameters  of  water  depth, 
wave  length  and  wave  height  were  varied  and  their  effect  on  forces 
determined.  The  horizontal  and  vertical  forces  on  the  model  were 
measured  by  instrumented  load  cells.  They  were  recorded  on 
electronic  recorders.  The  results  were  analyzed  with  the  aid  of  a 
digital  computer.  Forces  computed  using  existing  theories  were 
compared  to  measured  forces  and  the  results  used  to  determine 
inertial  coefficient.  The  results  of  the  data  were  plotted  in 
dimensionless  form  to  provide  a  correlation  between  ratios  of  wave 
length/water  depth,  wave  height/wave  length,  and  dimensionless 
force.  The  dimensionless  force  curves  can  be  used  to  estimate 
forces  due  to  waves  on  a  prototype. 
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INTRODUCTION 

In  recent  years,  offshore  oil  exploration  and  production  have 
moved  to  remote  areas  of  the  world  far  from  natural  or  artificial 
harbors.  Sometimes  the  harbors  in  the  remote  areas  do  not  have 
adequate  handling  capacities  to  accommodate  even  small  tankers. 
It  appears  that  savings  in  both  capital  investment  and  operating 
costs  may  be  realized  if  the  oil  could  be  stored  in  submerged 
tanks  in  the  oil  producing  areas.  This  conclusion  is  reached 
by  Chamberlin  (2)  in  his  consideration  of  underwater  tank  designs. 
Also  in  cases  of  small,  undeveloped  offshore  fields  it  may  not 
be  economical  to  pump  oil  usinq  a  pipeline  to  shore  and  to  develop 
shore  facilities. 

One  of  such  underwater  storage  tanks  was  recently  installed  in 
the  Middle  East  by  the  Chicago  Bridge  and  Iron  Company.  This 
tank  has  a  500,000  barrel  capacity,  is  about  260  feet  in  diameter 
and  about  65  feet  high.  It  was  installed  in  water  156-feet  deep. 
Since  such  a  tank  involves  a  large  investment  of  money  and  since 
no  adequate  theoretical  solutions  are  available  to  date,  extensive 
model  wave  tank  investigations  were  conducted  prior  to  design 
and  construction  of  the  tank  by  the  Chicago  Bridge  and  Iron 
Company. 

The  Institute  of  Ocean  Development  (8)  in  Tokyo  conducted  field 
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studies  in  Chiba  Port  adjacent  to  Tokyo  with  a  model  one-tenth 
the  size  of  the  prototype  tank.  The  upper  half  of  the  model  was 
hemispherical  in  shape  and  utilized  flexible  material.  The  lower 
half  of  the  model  was  elipsoidal  in  shape  and  was  made  from  steel. 
The  model  had  a  capacity  of  about  10,000  aallons. 

The  United  States  Commission  on  Marine  Science,  Engineering 
and  Resources  (3,  4)  and  the  Oceanographer  of  the  Navy  (14)  have 
predicted  that  underwater  facilities  will  be  built  at  an  increasing 
rate  in  the  near  future.  These  will  be  built  either  in  connection 
with  oil  and  gas  production,  or  in  connection  with  mineral  recovery 
from  the  ocean  bottom,  or  as  underwater  habitats.  Submerged  mili- 
tary support  bases  for  submersibles  are  predicted.  These  support 
bases  will  eliminate  sea  state  as  a  cause  for  the  prevention  or 
termination  of  operations  to  support  submersibles.  Habitats  can 
support  scientific  explorations  and  the  exploitation  of  the  natural 
resources  of  the  ocean  floor.  A  great  need  exists  to  be  able 
to  predict  forces  caused  by  waves  on  submerged  structures.  Until 
better  theoretical  force  predictions  are  available,  design  criteria 
will  have  to  rely  upon  model  tests. 

On  the  basis  of  model  tests,  this  thesis  provides  design  in- 
formation about  wave  forces  on  a  selected  number  of  shapes  of 
underwater  structures. 


LITERATURE  SURVEY 

Forces  on  a  body  oscillating  in  a  viscous  fluid  or  on  a  fixed 
body  in  an  oscillating  viscous  fluid  are  due  to  two  causes.  These 
are  inertia  and  drag.  The  drag  term  is  a  function  of  velocity, 
the  inertial  term  a  function  of  acceleration.  This  observation 
was  first  made  by  G.  G.  Stokes  (9)  in  1851. 

Forces  on  piles  due  to  gravity  waves  were  first  studied  by 
Morison,  O'Brien,  Johnson  and  Schaaf  (1)  in  1950.  They  set  forth 
the  so-called,  and  well  known  "Morison  Equation"  which  relates 
forces  to  wave  parameters.  The  force  was  found  to  consist  of 
two  components.  These  were  a  drag  force  proportional  to  the  square 
of  the  velocity  which  may  be  represented  by  a  drag  coefficient 
having  substantially  the  same  value  as  for  steady  flow  and  a 
virtual  mass  force  proportional  to  the  horizontal  component  of  the 
acceleration  force  exerted  on  the  mass  of  displaced  water.  The 
Morison  Equation  is: 


where 


9u 
=  rD  T  rM  =  l/a,D  p  M|U|U  T  UM  p  v 

1 


F  =  FD  +  FM  =  1/2CD  p  A|u|u  +  cM  p  V  |f (1) 


F  =  Total  force 
Fp  =  Drag  force 


The  literature  survey  will  use  the  symbols  as  they  appear  in 
the  cited  literature,  and  will  be  redefined  for  each  author.  Sym- 
bols for  other  sections  of  the  thesis  will  be  consistent  and  will 
be  defined  only  the  first  time  they  appear.  A  list  of  symbols 
is  qi ven  in  Appendix  I . 


FM  =  Inertial  force 
M 

Cn  =  Drag  coefficient 

CM  =  Inertial  coefficient 
M 

A  =  Cross  sectional  area  of  object  in  direction  perpen- 
dicular to  the  wave  travel 
V  =  Volume  of  object 
p  =  Mass  density 
u  =  Horizontal  velocity  of  water  particle 

r)U 

—r  =  Horizontal   acceleration  of  water  particle 

ot 

To  verify  the  theory  experimentally,  a  pile  was  attached  by  its 
base  to  the  bottom  of  a  wave  channel  and  restrained  by  a  load  cell 
at  the  top.  The  laboratory  results  showed  close  agreement  between 
computed  and  measured  force  with  the  exception  of  some  experimental 
scatter.  Further  work  on  waves  on  piles  was  published  by  Mori  son 
(11)  in  1951  and  Morison,  Johnson  and  O'Brien  (12)  in  1954 
utilizing  the  same  basic  relationship. 

In  1952  O'Brien  and  Morison  (13)  applied  the  Morison  Equation 
to  a  submerged  spherical  object.  The  components  of  force  were 
due  to  both  drag  and  inertia.  Forces  on  the  sphere  were  evaluated 
when  it  was  suspended  from  a  rod  above  the  bottom  and  when  it  was 
placed  on  the  bottom  by  use  of  a  tee  made  of  razor  blades.  The 
tee  was  made  by  attaching  two  razor  blades  to  the  tank  bottom. 
The  blades  were  perpendicular  to  the  direction  of  wave  travel  and 
parallel  to  each  other.  The  suspended  sphere  was  used  to  measure 
horizontal  force.  The  sphere  on  the  tee  was  used  to  measure 


vertical  force.  In  the  latter  case,  the  period  was  held  constant 
and  the  wave  height  raised  in  small  increments  until  the  sphere 
tipped  off  the  tee.  The  uplift  force  was  equal  to  the  weight  of 
the  sphere  when  movement  occured.  Drag  force  and  inertial  force 
were  computed  by  utilizing  the  measured  force  in  the  cases  where 
the  crest,  trough  and  nodes  of  the  wave  were  over  the  SDhere. 
In  1953  Reid  and  Bretschneider  (15)  utilized  the  Mori  son 
Equation  for  the  computation  of  wave  forces  on  piles  and  derived 
formulas  for  computation  of  wave  forces  on  large  submerged  objects 
The  forces  on  large  objects  were  found  to  be  almost  entirely 
inertial.  This  permitted  force  due  to  drag  to  be  neglected  in 
force  computations.  Due  to  the  fact  that  inertial  force  is  a 
function  of  particle  acceleration,  once  drag  is  neglected  the 
maximum  horizontal  force  occurs  when  the  wave  node  is  over  the 
model  center.  The  inertial  force  was  computed  directly  from  the 
horizontal  distribution  of  pressure.  The  equations  developed  by 
Reid  are  simple,  straight-forward  and  accurate  for  large  objects 
where  drag  forces  are  negligible  and  inertial  forces  predominant. 
For  forces  on  a  large  submerged  object  utilizinq  small  amplitude 
(Airy)  wave  theory: 

FH  =  CML2L3  M  [cose.,  -  cos  (e1  +  -—-)] (2) 


where: 


Fn  =  Total   horizontal    force 
CM  =   Inertial   coefficient 


L,  =  Length  of  structure  (in  direction  of  wave  travel) 

L?  =  Width  of  structure 

U  =  Height  of  structure 

K  =  Pressure  factor  = 


cosh  2-ttcI 
L 
d  =  Water  depth 

L  =  Wave  length 

y     =  Specific  weight 

H     =  Wave  height 

e,  =  Phase  angle  (position)  of  end  of  model  facing  waves 

measured  from  crest  in  direction  of  wave  travel 

The  maximum  value  of  F  for  this  wave  will  occur  when  the 

center  of  the  object  is  at  the  position  e  =  -S-  radians.  Thus: 

FMax  "  CHL2L3  K"   H  s1n  ^ (3) 

A  re-analysis  of  all  existing  data  on  wave  forces  on  piles 
was  made  by  Crooke  (5)  in  1955.  His  paper  reconciled  differences 
in  previous  studies  by  applying  a  consistent  approach  which 
resulted  in  reasonable  results.  He  applied  the  Iverson  approach 
to  wave  forces.  The  analysis  was  made  utilizinq  only  one  co- 
efficient in  the  wave  force  equation.  This  coefficient  was 
dependent  both  upon  velocity  and  acceleration. 

ai   VL     \/2 

C  =  <p     (%,  -*-,     U (4) 

yd.        v  gl_ 


where: 

AL 

— 7r      =  Iversen's  modulus 

T 

VL 

— -    =  Reynolds  number 

-t-r      =  Froude  number 

gL 

A    =  Particle  acceleration 

V    =  Particle  velocity 

L    =  Characteristic  length 

p    =  Density 

v         =   Viscosity 

g    =  Gravity 
Crooke's  approach  worked  well  for  deep  and  shallow  water  waves 
as  long  as  the  appropriate  wave  theory  was  used. 

In  1958  Keulegan  and  Carpenter  (9)  investigated  the  inertial 
and  drag  coefficients  of  cylinders  and  plates  when  exposed  to 
gravity  waves.  They  found  that  a  remainder  function  AR  could 
be  introduced  to  provide  a  truer  representation  of  force  when  con- 
sidering the  coefficients  CM  and  CD  as  being  constant  through 
the  wave  cycle.  AR  was  required  because  the  point  values  of 
CM  and  Cn  differed  from  their  average  values.  They  also  found 

that  average  values  of  Cn  and  CM  could  be  correlated  with  point 

UmT 
values  in  the  MoriSon  Equation  by  introducing  a  parameter  -sp  . 


where: 


U     =  Maximum  horizontal   particle  velocity 
m 


T  =  Wave  period 

D  =  Diameter  of  cylinder  or  width  of  plate 
The  Mori  son  Equation  thus  modified  becomes: 


F 
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p  U  mD  m 

AR  is  obtained  by  subtracting  the  computed  values  for 

r 

/TT   r   x   /Do     .x        D     1        i 

(4  CM)  (y-sine)  -  -j     |cose|cose 
m 

from  observed  values  for 


*UmD 
It  was  found  that  inertial  coefficient  and  drag  coefficient  varied 
with  the  diameter  of  the  cylinder.  For  a  given  cylinder  it 

varied  with  the  values  of  maximum  current  U  .  A  strong  correla- 

U  T     m 
tion  between  the  period  parameter  -^  and  inertial  and  drag  coef- 


ficients  was  found  to  exist. 

It  was  found  that  for  the  plates  the  hiqher  values  of  the  drag 
coefficient  are  obtained  at  the  smaller  values  of 

UmT 

m 
D 

and  the  higher  values  of  inertia!  coefficient  are  associated  with 

larger  values  of 

U  T 
m 

D 
Brater,  McNown  and  Stair  (1)  investigated  experimentally  the 
effect  of  gravity  waves  on  models  of  submerged  barge-like  structures 
in  1958.  They  evaluated  a  flat  plate,  a  right  parallelepiped  and 
two  parallelepipeds  with  a  somewhat  modified  shape.  The  models 
were  suspended  from  rods  instrumented  with  load  cells.  The 
Morison  Equation  was  utilized  for  analysis. 

F  =  1/2  CD  A  p  |  u  |  u  +  Cm  p  V  |£ (7) 


where 


Cn  =  Drag  coefficient 

A  =  Cross  sectional  area  of  model 

p  =  Density 

u  =  Horizontal  particle  velocity 


C  =  Inertial  coefficient 
m 


V  =  Volume  of  model 


9U 

rr  =  Horizontal  acceleration  of  narticle 


10 


x 

The  maximum  force  on  the  model  structures  occurred  between  p  =  0 
and  £  =  1/4 

where: 

L  =  Wave  length 

X  =  Distance  from  wave  crest  to  location  of  maximum 
force  measured  in  direction  of  wave  travel 

Drag  forces  were  found  to  be  negligible  on  the  barge.  Analysis 
of  the  results  were  made  utilizing  the  inertial  portion  only.  The 
analysis  of  inertial  force  used  pressure  to  evaluate  horizontal 
force. 

<P1  -  V  \  ■  >  V  ft <8> 

where: 

P-,  =  Pressure  on  front  of  model 
P?  =  Pressure  on  rear  of  model 

This  expression  is  substituted  in  the  inertial  portion  of  the 
Morison  Equation  to  yield 


Fi  -  Cm  <P1  "  P2>  Ax <9> 


where: 


F.  =  Portion  of  force  due  to  inertia 


The  equation  of  motion  in  the  vertical  direction  is 

dp                           9V  /,nv 

DZ  =  "W  +  pJt <10) 


11 


where: 

t-2-  =   Increment  of  pressure 

9Z  K 

w     =  Specific  weight  of  water 
When  Eq.    10  is   substituted  into  the  Airy  expressions   for  particle 
acceleration  and  integrated,   the  result  is 

p  =  W(.z  +  y,  +  0  I!"     cosh  Md+z)/L  -^cosh  2„(d+y)/L  cos6   (n) 

where: 

y  =  Elevation  of  surface  above  still  water 

z  =  Distance  from  water  surface  to  model 

H  =  Wave  height 

T  =  Wave  period 

d  =  Depth  of  water 
For  the  case  of  maximum  force  which  occurs  when  the  model  is 
centered  at  the  quarter  point  of  the  wave,  the  above  equation  be- 
comes : 

2ttX,  9      ttH    cos    2ttX,/L 

P1   -  P2  =  wH  cos  ~-  [K  -   2  sinli     ( ^ ^    ■    ■      (12) 

where: 

x,  =  Horizontal  distance  from  wave  crest  to  leadinq 

edge  of  barge 

K  =  cosh  27r(d+z)/L 
cosh  2ttc!/L 

When  measured  values  are  substituted  in  Eq.  12  the  value  of 
(P,  -  P?)  can  be  computed.  Eq .  9  can  be  rewritten  as 

Cm  =  (P]  -  P2)Ax (13) 


12 


from  which  C  can  be  computed. 

m 

Both  drag  force  and  inertial  force  were  found  to  be  significant 
in  the  results  from  the  flat  plate.  Values  of  Cn  and  CM  which 
yielded  the  best  correspondence  with  measured  forces  and  phase 
angles  were  utilized  in  the  Morison  Equation.  The  values  for 
computed  wave  parameters  were  derived  utilizing  Airy  theory.  For 
any  given  wave,  the  force  is  a  function  of  the  phase  angle  only. 
The  Morison  Equation  was  modified  as  follows  to  derive  the  phase 
angle  for  maximum  force  where  drag  and  inertial  forces  were  both 
significant, 

F  =  C  '  Cos2  e  +  C  '  Sin  e (14) 

Dm  v 

where: 

F  =  Total  force 

e  =  Phase  anqle 


c  •  =  r  A  p  ^4  [C°Sh  2\{i+Z)Jh2 (15) 

D  D  2T2     L         sinh  2u  d/LJ 

r    >    -  r  \t  2tt  H     r-cosh  2tt   (d+z)/Ln  /-,,-\ 

Cm     "  Cm  p  V  -^2—  [     sinh  2,  d/L     ] (16) 

The  phase  angle  for  maximum  force  is  obtained  by  setting  dF/de  =  0 
with  the  result 

si"W^ <17> 

They  found  that  inertial  coefficient  for  the  barge  was  1.5. 


13 


Forces  obtained  from  tests  of  a  barge  with  rounded  edaes  and  a 

barge  with  a  slot  in  the  center  were  not  substantially  different 

from  those  obtained  from  tests  on  a  normal  barge  with  square 

edges.  This  was  due  to  the  fact  that  inertial  force  is  a  function 

of  volume.  Measured  inertial  forces  on  the  barge  agrees  well  with 

theory.  The  drag  coefficient  for  the  plate  was  3.5  and  virtual 

mass  (the  product  of  volume  and  C  )  was  1.75  times  that  of  a 

m 

circumscribing  cylinder. 

A  recent  paper  by  Evans  (6)  published  results  of  study  of 
forces  on  piles  of  full -size  prototype  offshore  platforms.  He 
utilized  the  Mori  son  Equation  to  compare  computed  predictions  with 
measured  results.  The  measured  data  included  waves  of  large  size 
(up  to  106  ft.)  due  to  hurricanes.  His  data  was  taken  during  two 
extended  periods.  They  were  1954  to  1958  and  1960  to  1963.  The 
study  proved  the  accuracy  of  the  Mori  son  Equation.  The  average 
values  of  measured  and  theoretical  forces  were  within  15  percent 
of  each  other.  The  value  obtained  for  Cn  was  0.5.  The  value 
for  C  was  1.5.  This  agrees  well  with  previous  work  done  by 
others  and  summarized  by  Wilson  and  Reid  (17).  The  drag  coef- 
ficients were  found  to  decrease  for  increasina  wave  height.  The 
inertial  coefficient  did  not  change  greatly  with  moderate  changes 
in  wave  size.  However,  it  doubled  in  value  from  a  23  ft.  wave  to 
a  106  ft.  wave. 
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THFORFTICAL  CDNSIDFRATIOMS 

Progressive  linear  wave  (Airy)  theory  is  used  for  this  study. 
A  brief  presentation  is  made  of  the  portions  aDproDriate  to  this 
work.  Detailed  development  is  made  by  Faaleson  and  Dean  (7)  and 
Wiegel  (16). 

Airy  theory  makes  the  followina  assumptions: 

That  the  bottom  is  smooth  and  impermeable 

That  the  fluid  is  incompressible,  irrotational  and  inviscid 

That  the  depth  is  constant  (measured  from  water  level) 

That  surface  tension  is  nealected 

That  cori olis  force  is  neglected 

That  atmospheric  pressure  is  constant 

That  the  height  of  the  surface  is  small  in  comparison  with 
wave  lenath 

That  all  non-linear  terms  may  be  neglected 

The  wave  surface  is  described  by  the  following  relationship: 

n  =  \    cos  2tt  (£  -  |) (1R) 

where: 

n  =  Distance  of  the  free  surface  from  mean  water  level 

x  =  Distance  from  the  wave  crest  measured  in  the  direction 

of  wave  travel 
L  =  Wave  length 
H  =  Wave  height 
t  =  Time 
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T  =  Wave  period 

Some  of  the  various  wave  measurements  and  parameters  which  will 
be  used  in  this  section  and  the  results  and  conclusions  section 
are  shown  in  Fig.    1 . 

The  velocity  potential   for  the  small   amplitude  wave  is  given 

by  the  expression 

_  ga     cosh  k   (z+d)     „rtc    ,.  .%  nQ>> 

<i>  =  •*!—    1    ■  \ *-    cos   (kx  -  at) { iy; 

^         a       cosh  kd  v  ' 

where: 

4>  =  Velocity  potential 

g  =  Gravity 

a  =  Wave  amplitude  =  H/2 

o  =  2tt/T 

z  =  Distance  from  still  water  level 

d  =  Water  depth 

k  =  2tt/L 

The  celerity  of  the  wave  is  given  by  the  expression: 

C  =  51  tanh  %±    =  [Sk  tanh  ^]1/2 (20) 

where: 

C  =  Celerity 
The  expression  for  wave  length  is: 

L  =  5li  tanh  M (21) 

In  deep  water  the  depth  becomes  large  relative  to  wave  length 
and  the  hyperbolic  tangent  function  (tanh  ——)   approaches  unity. 
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Wave  celerity  becomes 

C0  ■  [f^]1/2 (22) 

where: 

C     =  Deep  water  celerity  and  wave  length  becomes 

nT2  7 

|_     =  2i_    =  5.12  r (23) 

O  tTT 

where: 

L  =  Deep  water  wave  length 

The  wave  period  is  an  independent  variable.  It  also  remains 
relatively  constant  irregardless  of  changes  in  wave  height,  water 
depth,  and  wave  length.  If  corresponding  expressions  for  general 
waves  and  deep  water  waves  are  solved  using  a  constant  depth,  the 
following  expression  results 

d  tanh  2^d  =  d_ (24) 

o 

This  enables  direct  solution  for  the  wave  length  from  the  deep 
water  wave  length  which  is  easily  calculated  from  Eq.  23. 

The  horizontal  component  of  water  particle  velocity  is  given 
by  the  expression. 

»•&-"£  £2a#frI  cos(kx-0t) (25) 

where: 

u  =  Horizontal  velocity  of  water  particle 
The  horizontal  component  of  water  particle  acceleration  is 
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given  by  the  expression: 

du           ,  cosh  k(z+d)   .   /,     .v  ,nc\ 

U  =  3t  =  gak   cosh  kd    sin  (kx  -  at) (26) 

where: 

u  =  Horizontal  acceleration  of  water  particle 
The  Morison  Equation  relates  force  to  particle  motion. 

Fh  =  1/2  p  CD  A|u|u  +  p  Cm  V  u (27) 

where: 

p  =  Density 

F,  =  Horizontal  force 

For  a  given  vertical  increment  of  the  model,  this  expression  be- 
comes 

dFh  =  dFD  +  dFj (28) 

where: 

FD  =  Force  due  to  drag 
Fj  =  Force  due  to  inertia 

dFT   =  C     p  WL     u  dz (29) 

I         m         m  m  v/ 

where: 

W  =  Width  of  model 
m 

L  =  Length  of  model 
m     3 

Hm  =  Weight  of  model 
m     3 

The  total  inertial  force  can  be  obtained  from  substitution  of 

Eq.  26,  in  Eq.  29  and  integration  of  the  resultant  expression. 
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FT   =  C     p   W  L 
I         m         m  m 


-d  +  H 


m 


2tt  h  cosh  k(d+z)      .„      .     .  ,,n* 

"  ~p     sinh  kd sin  at  dz-    '(30) 


(cmpwmLm)c2; H  sin  at) 

m  m  m   T2  sinh  kd 


-d  +  H 


m 


cosh  k(d+z)  dz 


fr     .,  ,    \/-2tt  H  sin  atwsinh  k(d+z) 

(CmpWmLm)  (  T2  ginh  -H   }  ( IT L 

sinh  kd 


-d  +  H 


m 


-d 


?.2    H    <Hn    „t       Sinh    k(d-d+Hm)    "    Sinh    k(d"d) 

_    /r     u  i     w^-71     "i  sin  otw  m 

-   (C  PW  L    )(— ^ )(— £ 

m    m  m      r  sinh  kd  K 


9      2U      ef«  +  Sinh      dk      Hm 

/r     ii  i     \c^  rl  sin  atw                   rru 
(C  PW  L    )  (— )  ( r ) 

m  m  m   r  sinh  kd      K 


(31) 


The  force  due  to  drag  on  a  vertical  increment  of  model  is 


dFn  =  1/2  Cn  p  Wju|u  dz 


D 


D  w   "m' 


(32) 


The  total  drag  force  can  be  obtained  by  substitution  of  Eq.  25  in 
Eq.  32  and  integration  of  the  resultant  expression. 

-d  +  H 


m 


FD  *  V2  CD  "  \ 


ulu  dz 


(33) 
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1/2  CD  p  Wm 


■d  +  H 


m 


/ttHx2     ,cosh  k(d+z)   „rte      .  ,2    . 
(-T}        (     sinh  kd       cos  ot)     dz 


,  /0  r          ..      /ttH\2      /  I  cos   otjcos   at\ 
1/2    Cn    p    W       (-y)  (J L^-* ) 

D         m       '  (sinh  kdT 


■d  +  H 


m 


cosh     k(d+z)  dz 


,  /0  r         ,,     /ttHx2   /Icos  atlcos  at\rsinh  2k(d+z)       k(d+z)n 

1/2  Cn  P  W  (-y)  (J 1 s — )L r —  +    ?  ■* 

D    m   T     (sinh  kdT 


-d  +  H 


-d 


m 


1/2  r  P  W  (^)2  (Icos  at|cos  at)( 
D    m   T      (sinh  kdT 


sinh  2k(d-d+H  )   k(d-d+H  ) 


sinh  2k(d-d)    k(d-d)  * 
4      ~         2       ' 


u  o  i     .  i     +  sinh  2k  H    k  H 

t  ,o  r         H  fT*n\<-   /Cos  ot  cos  ot\t  m  ,  ni\ 

1/2  cn  p  wm  ^     < ? — ^ 3 — T~> 

D   m   '  (sinh  kdT       a       l 


.  (34) 


If  a  submerged  object  has  a  large  volume,  the  inertial  forces 
predominate,  and  the  drag  forces  can  be  ignored.  If  it  has  a 
length  in  the  direction  of  wave  travel  which  is  large  compared  to 
the  wave  length,  an  analysis  of  force  can  be  computed  directly 
from  the  horizontal  distribution  of  pressure  beneath  the  wave. 
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The  difference  between  the  actual  gauge  pressure  at  the  bottom 
and  the  calculated  hydrostatic  pressure  based  on  the  still  water 
depth  is  defined  as  pressure  anomaly. 

AP  =  P  -  Yd (35) 

where: 

aP  =  Pressure  anomaly 

P  =  Gauge  pressure 

Y  =  Specific  weight  of  water 
A  pressure  factor  which  is  a  function  of  H/d  is  defined  as  follows 

Up)  -  UpL 
K  =  *fl " <36> 

where: 

K  =  Pressure  factor 

Up)  =  Pressure  anomaly  under  the  wave  crest 

Up).  =  Pressure  anomaly  under  the  wave  trough 

k  = Vrr (37) 

cosh  (-j— ) 

For  Airy  wave  theory  the  variation  of  Ap/yH  with  phase  angle  is  a 
function  (-^)  cos  e  where  e  is  the  phase  angle,  and  is  equal  to 
ot  or  kx. 

The  general  equation  for  horizontal  force  acting  on  a  large 
rectangular  submerged  object  is 

F.  =  C  H  W  [AP,  -  AP0] (38) 

h    m  m  m        2J  v  ' 
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where: 

aP,  =  Pressure  anomaly  at  leading  edge  of  object 
APp  =  Pressure  anomaly  at  trailing  edge  of  object 

In  the  special  case  of  the  Airy  theory  wave 

Fh     =     Cm     Hm    Wm       ^T        CCOS01      "     C0S(61     +     2tT    7^)1 (39) 

n        m    m    m      2  1  1  L 

where: 

9l   =  2ttX1 


L 

x,  =  Horizontal  distance  from  the  leadina  edge  in  direc- 
tion of  wave  travel  relative  to  a  crest 

The  maximum  value  of  the  horizontal  force  (Fu)    occurs  when  the 

H  max 

center  of  the  tank  is  at  the  position 

6       2 
radians.     Therefore, 

<FH>max  *  Cm  Hm  Wm  ^  H  si"   hr> <40> 
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EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

Models 

With  the  exception  of  the  flat  plate  the  models  were  all 
constructed  from  plexiglass.  The  flat  plate  was  fabricated 
from  aluminum.  All  models  used  (with  the  exception  of  the  aluminum 
plate)  are  shown  in  Fig.  2  and  Fig.  3.  The  models  consisted  of 

(a)  plate  7-23/32  in.  wide,  4-1/2  in.  high  and  5/8  in.  thick. 

(b)  rectangular  tank  7-23/32  in.  wide,  4-1/2  in.  high  and  2  in. 
long. 

(c)  rectangular  tank  7-23/32  in.  wide,  4-1/2  in.  high  and  4 
in.  long. 

(d)  rectangular  tank  7-23/32  in.  wide,  4-1/2  in.  high  and 
8  in.  long. 

(e)  a  half-cylinder  tank  having  radius  equal  to  4-15/32  in. 
and  length  equal  to  7-11/16  in. 

Equipment  and  Procedure 

The  models  were  placed  in  a  three  feet  deep,  two  feet  wide  and 
one  hundred  and  twenty  feet  long  wave  channel.  Waves  were  generated 
with  an  oscillating  pendulum- type  generator  capable  of  producing 
a  wide  range  of  wave  lengths  and  heights.  The  generator  is  powered 
with  a  variable  speed  motor.  The  wave  period  of  the  waves  was 
varied  by  changing  the  speed  of  the  motor.  The  wave  generator 
is  shown  in  Fig.  4. 

The  wave  height  was  changed  by  varying  the  adjustment  on  the 
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FIG.  2.-MnnFLS  STUDIED 


FIG.  3. -CLOSE-UP  OF  FOUR  INCH  LONG  RECTANGULAR  MODEL 
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FIG.  4. -WAVE  GENERATOR 


FIG.  5. -EXPERIMENTAL  FACILITY 
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excentricity  of  the  paddle  arm  on  the  flywheel  which  resulted  in 
a  change  in  the  stroke  of  the  pendulum.  The  pendulum  was  adjusted 
to  provide  the  proper  shape  for  each  wave.  A  stroke  which  had 
equal  amplitude  from  top  to  bottom  was  used  for  the  most  shallow 
water  waves.  A  stroke  with  a  large  amplitude  at  the  top  and  a  small 
amplitude  at  the  bottom  was  used  for  deep  water  waves.  An  efficient 
absorber  is  installed  at  the  other  extremity  of  the  tank  to  prevent 
reflected  waves  from  interfering  with  generated  incident  waves. 
The  absorber  has  a  permeable  face  with  an  impermeable  backina.  The 
channel  has  3/8  inch  glass  side  panels  throughout  its  length  to 
allow  for  visual  observation  of  experiments. 

The  experiments  were  conducted  at  water  depths  of  24  inches, 
18  inches  and  13  inches.  Wave  lengths  were  varied  to  produce 
several  L/d  ratios  between  1.5  and  6.0.  This  covers  deep  water 
waves  and  the  deep  half  of  intermediate  waves  (where  intermediate 
waves  are  defined  as  having  a  L/d  ratio  between  2  and  20).  The 
pendulum  was  timed  by  stopwatch  through  20  to  50  cycles  and  the 
period  determined  to  the  closest  1/100  second. 

Approximately  100  different  combinations  of  water  depth,  wave 
period  (and  consequently  wave  length)  and  wave  heiqht  were  used  for 
each  model.  The  wave  heights  were  increased  until  the  waves  broke. 
The  schedule  of  test  waves  used  are  shown  in  Table  1.  For  each 
combination  of  water  depth  and  wave  period  the  schedule  increases 
wave  height  to  a  value  which  will  yield  a  steepness  equal  to  the 
theoretical  maximum  of  1.43.  A  total  of  599  different  experimental 
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TABLE   1 


TEST  WAVE  SCHEDULE 


Test 

Depth 

Period 

Wave 

d/L 

L/d 

H/L 

Wave 

in 

in 

Height 

Number 

feet 

seconds 

in   feet 

1 

2.000 

1.68 

0.10 

0.178 

5.61 

0.009 

2 

2.000 

1.68 

0.20 

0.178 

5.61 

0.018 

3 

2.000 

1.68 

0.30 

0.178 

5.61 

0.027 

4 

2.000 

1.68 

0.40 

0.178 

5.61 

0.036 

5 

2.000 

1.68 

0.50 

0.178 

5.61 

0.044 

6 

2.000 

1  .68 

0.60 

0.178 

5.61 

0.053 

7 

2.000 

1.68 

0.70 

0.178 

5.61 

0.062 

8 

2.000 

1.55 

0.10 

0.203 

4.93 

0.010 

9 

2.000 

1.55 

0.20 

0.203 

4.93 

0.020 

10 

2.000 

1.55 

0.30 

0.203 

4.93 

0.030 

11 

2.000 

1.55 

0.40 

0.203 

4.93 

0.041 

12 

2.000 

1.55 

0.50 

0.203 

4.93 

0.051 

13 

2.000 

1.55 

0.60 

0.203 

4.93 

0.061 

14 

2.000 

1.55 

0.70 

0.203 

4.93 

0.071 

15 

2.000 

1.29 

0.10 

0.265 

3.77 

0.013 

16 

2.000 

1.29 

0.20 

0.265 

3.77 

0.027 

17 

2.000 

1.29 

0.30 

0.265 

3.77 

0.040 

18 

2.000 

1.29 

0.40 

0.265 

3.77 

0.053 

19 

2.000 

1.29 

0.50 

0.265 

3.77 

0.066 

20 

2.000 

1.29 

0.60 

0.265 

3.77 

0.080 

21 

2.000 

1.29 

0.70 

0.265 

3.77 

0.092 

22 

2.000 

1.09 

0.10 

0.342 

2.92 

0.017 

23 

2.000 

1.09 

0.20 

0.342 

2.92 

0.034 

24 

2.000 

1.09 

0.30 

0.342 

2.92 

0.051 

25 

2.000 

1.09 

0.40 

0.342 

2.92 

0.068 

26 

2.000 

1.09 

0.50 

0.342 

2.92 

0.086 

27 

2.000 

1.09 

0.60 

0.342 

2.92 

0.103 

28 

2.000 

0.89 

0.10 

0.505 

1.98 

0.025 

29 

2.000 

0.89 

0.20 

0.505 

1.98 

0.051 

30 

2.000 

0.89 

0.30 

0.505 

1.98 

0.076 

31 

2.000 

0.89 

0.40 

0.505 

1.98 

0.101 

32 

2.000 

0.89 

0.50 

0.505 

1.98 

0.126 
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TABLE    1 .    - 

TEST  WAVE 

SCHEDULE 

(Continued) 

Test 

Depth 

Period 

Wave 

d/L 

L/d 

H/L 

Wave 

in 

in 

Height 

Number 

feet 

seconds 

in  feet 

33 

2.000 

0.78 

0.10 

0.646 

1  .55 

0.032 

34 

2.000 

0.78 

0.20 

0.646 

1.55 

0.065 

35 

2.000 

0.78 

0.30 

0.646 

1.55 

0.097 

36 

2.000 

0.78 

0.40 

0.646 

1.55 

0.129 

37 

2.000 

0.78 

0.50 

0.646 

1.55 

0.161 

38 

1.500 

0.71 

0.10 

0.579 

1.73 

0.039 

39 

1.500 

0.71 

0.20 

0.579 

1.73 

0.077 

40 

1.500 

0.71 

0.30 

0.579 

1.73 

0.116 

41 

1.500 

0.71 

0.40 

0.579 

1.73 

0.154 

42 

1.500 

0.79 

0.10 

0.467 

2.14 

0.031 

43 

1.500 

0.79 

0.20 

0.467 

2.14 

0.062 

44 

1.500 

0.79 

0.30 

0.467 

2.14 

0.093 

45 

1.500 

0.79 

0.40 

0.467 

2.14 

0.125 

46 

1.500 

0.79 

0.50 

0.467 

2.14 

0.156 

47 

1.500 

0.92 

0.10 

0.353 

2.84 

0.023 

48 

1.500 

0.92 

0.20 

0.353 

2.84 

0.047 

49 

1.500 

0.92 

0.30 

0.353 

2.84 

0.070 

50 

1.500 

0.92 

0.40 

0.353 

2.84 

0.094 

51 

1.500 

0.92 

0.50 

0.353 

2.84 

0.117 

52 

1.500 

0.92 

0.60 

0.353 

2.84 

0.140 

53 

1  .500 

1.09 

0.10 

0.265 

3.78 

0.018 

54 

1.500 

1.09 

0.20 

0.265 

3.78 

0.035 

55 

1.500 

1.09 

0.30 

0.265 

3.78 

0.053 

56 

1.500 

1.09 

0.40 

0.265 

3.78 

0.071 

57 

1.500 

1.09 

0.50 

0.265 

3.78 

0.088 

58 

1.500 

1.09 

0.60 

0.265 

3.78 

0.106 

59 

1.500 

1.31 

0.10 

0.200 

5.01 

0.013 

60 

1.500 

1.31 

0.20 

0.200 

5.01 

0.027 

61 

1.500 

1.31 

0.30 

0.200 

5.01 

0.040 

62 

1.500 

1.31 

0.40 

0.200 

5.01 

0.053 

63 

1.500 

1.31 

0.50 

0.200 

5.01 

0.067 

64 

1.500 

1.31 

0.60 

0.200 

5.01 

0.080 

65 

1.500 

1.31 

0.70 

0.200 

5.01 

0.093 

TABLE   1.    -  TEST   WAVE   SCHEDULE 
(Continued) 
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Test 

Depth 

Period 

Wave 

d/L 

L/d 

H/L 

Wave 

in 

in 

Height 

Number 

feet 

seconds 

in  feet 

66 

1.500 

1.47 

0.10 

0.171 

5.85 

0.011 

67 

1.500 

1.47 

0.20 

0.171 

5.85 

0.023 

68 

1.500 

1.47 

0.30 

0.171 

5.85 

0.034 

69 

1.500 

1.47 

0.40 

0.171 

5.85 

0.046 

70 

1.500 

1.47 

0.50 

0.171 

5.85 

0.060 

71 

1.500 

1.47 

0.60 

0.171 

5.85 

0.068 

72 

1.500 

1.47 

0.70 

0.171 

5.85 

0.080 

73 

1.083 

1.26 

0.10 

0.169 

5.93 

0.016 

74 

1.083 

1.26 

0.20 

0.169 

5.93 

0.031 

75 

1.083 

1.26 

0.30 

0.169 

5.93 

0.047 

76 

1.083 

1.26 

0.40 

0.169 

5.93 

0.062 

77 

1.083 

1.26 

0.50 

0.169 

5.93 

0.078 

78 

1.083 

1.26 

0.60 

0.169 

5.93 

0.093 

79 

1.083 

1.04 

0.10 

0.220 

4.54 

0.020 

80 

1.083 

1.04 

0.20 

0.220 

4.54 

0.041 

81 

1.083 

1.04 

0.30 

0.220 

4.54 

0.061 

82 

1.083 

1.04 

0.40 

0.220 

4.54 

0.081 

83 

1.083 

1.04 

0.50 

0.220 

4.54 

0.102 

84 

1.083 

1.04 

0.60 

0.220 

4.54 

0.122 

85 

1.083 

0.89 

0.10 

0.283 

3.53 

0.026 

86 

1.083 

0.89 

0.20 

0.283 

3.53 

0.052 

87 

1.083 

0.89 

0.30 

0.283 

3.53 

0.078 

88 

1  .083 

0.89 

0.40 

0.283 

3.53 

0.105 

89 

1  .083 

0.89 

0.50 

0.283 

3.53 

0.131 

90 

1.083 

0.89 

0.60 

0.283 

3.53 

0.157 

91 

1.083 

0.75 

0.10 

0.382 

2.62 

0.035 

92 

1.083 

0.75 

0.20 

0.382 

2.62 

0.070 

93 

1  .083 

0.75 

0.30 

0.382 

2.62 

0.106 

94 

1.083 

0.75 

0.40 

0.382 

2.62 

0.141 

95 

1.083 

0.71 

0.10 

0.421 

2.37 

0.039 

96 

1.083 

0.71 

0.20 

0.421 

2.37 

0.078 

97 

1.083 

0.71 

0.30 

0.421 

2.37 

0.117 

98 

1.083 

0.71 

0.40 

0.421 

2.37 

0.156 
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runs  were  made.  The  resulting  massive  amount  of  data  provided 
almost  total  coverage  within  the  range  of  wave  parameters  tested. 

An  over-all  view  of  the  experimental  facility  is  shown  in 
Fig.  5.  A  schematic  diagram  of  the  experimental  facility  is  shown 
in  Fig.  6.  To  accomplish  the  force  measurements  the  model  was 
supported  by  four  vertical  wires  with  the  upper  ends  of  the  wires 
supported  from  canti levered  load  cells.  The  wire  was  0.015  inch 
diameter  stainless  steel  cable.  It  was  composed  of  7  strands  and 
provided  a  supple,  inelastic  support  for  the  model.  It  was 
fastened  to  the  model  by  brass  and  stainless  steel  swivel  snap 
hooks  at  the  lower  end  and  connected  to  the  load  cell  by  a  threaded 
eye  bolt  at  the  upper  end.  The  cable  was  fastened  to  the  snap 
hook  and  threaded  eye  by  use  of  a  crimped  stainless  steel  sleeve. 
The  threaded  eye  had  1  1/2  inches  of  threaded  shaft  to  facilitate 
vertical  adjustment  of  the  models.  The  models  were  placed  with 
their  base  1/4  inch  above  the  tank  bottom.  A  close-up  view  of 
the  experimental  set  up  is  shown  in  Fig.  7. 

Vertical  forces  caused  the  load  cells  to  move  in  a  vertical 
direction.  The  load  was  measured  by  strain  gauges  and  recorded 
on  electronic  recorders. 

Since  wave  forces  produce  loads  both  ud  and  down  on  the  model, 
it  was  necessary  to  have  it  weighted  with  metal  inside  the  model 
so  that  the  vertical  load  lines  were  pretensioned.  This  pre- 
tension was  of  high  enough  value  so  that  the  lines  never  went 
slack.  By  hanging  the  structure  from  four  load  lines  it  was  also 
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FIG.    7. -CLOSE-UP  OF  EXPERIMENTAL  SET-UP 
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possible  to  measure  moments  on  the  model.  A  model  suspended  from 
the  vertical  load  cells  is  shown  in  Fig.  8. 

The  horizontal  forces  were  measured  as  follows:  A  wire  and 
spring  were  used  to  attach  the  forward  end  of  the  model  to  the  wave 
tank  bottom.  The  downstream  end  of  the  model  was  attached  by  a 
wire  which  was  diverted  ninety  degrees  over  a  pulley  to  a  load  cell. 
The  spring  provided  a  constant  tension,  the  variation  of  which  was 
measured  by  the  load  cell  and  indicated  the  wave  induced  horizontal 
force  on  the  model.  The  measured  load  was  recorded  on  an  electronic 
recorder. 

Both  the  stud  and  the  pulley  were  placed  at  a  horizontal  dis- 
tance of  three  feet  from  the  model  to  minimize  disturbance  in  the 
model  area.  A  large  diameter  (5  inches)  ball  bearing  pulley  was 
used  to  minimize  the  effects  of  friction.  A  schematic  diagram  of 
the  model  and  load  cells  is  shown  in  Fig.  9.  Known  weights  were 
used  to  calibrate  the  force  measurement  apparatus.  A  capacitance 
type  wave  gauge  was  used  to  measure  wave  profiles.  The  wave  gauge 
was  calibrated  by  immersing  the  gauge  wire  into  the  water  by  incre- 
mented known  distances  with  the  water  still  in  the  channel.  The 
wave  height  and  wave  profiles  were  recorded  on  an  electronic 
recorder. 

A  probe  was  positioned  adjacent  to  the  model  so  that  the 
passing  of  each  wave  crest  past  the  model  would  be  marked  on  the 
recorder.  This  made  possible  the  measurement  of  the  phase  angle 
(or  time  difference)  between  the  crest  and  the  maximum  wave  forces. 
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FIG.  8. -HALF-CYLINDER  MODEL  IN  PLACE 
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A  synchronized  marker  was  used  to  place  a  common  time  mark  simul- 
taneously on  each  recorder  data  paper.  This  enabled  the  data 
outputs  from  the  three  recorders  to  be  correlated  in  time.  The 
data  was  recorded  on  hot  needle- type  direct-recording  equipment. 
The  wave  height  and  wave  profiles  were  recorded  on  a  Hewlett-Packard 
dual -channel  carrier  pre-amplifier  recorder  (model  321).  The 
wave  forces  were  recorded  on  Sanborn  two-and  four-channel  recorders 
(model  150).  Recorder  outputs  are  shown  in  Fig.  10  and  Fig.  11. 
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PRESENTATION  AND  DISCUSSION  OF  RESULTS 

The  primary  objective  of  the  research  was  to  determine  ex- 
perimentally the  forces  due  to  gravity  waves  on  several  submerged 
structures  of  rectangular  and  half  cylindrical  shape.  The  relation 
between  measured  vertical  and  horizontal  force  and  various  wave 
parameters  was  investigated.  Forces  determined  by  experiment  were 
compared  to  those  calculated  using  existing  theories  for  wave 
forces  on  submerged  objects.  Wave  reflections  from  the  model  were 
assumed  to  be  small.  Therefore,  they  were  not  measured  and  were 
neglected  in  the  theoretical  analysis. 

The  raw  data  was  in  the  form  of  recorder  output  shown  in 
Fig.  10  and  Fig.  11.  The  period  was  determined  by  measurement 
of  the  length  in  decimeters  between  a  wave  crest,  and  the  crest 
of  the  20th  wave  following  it.  The  paper  feed  was  5  millimeters 
per  second,  so  twenty  waves  allowed  reading  the  measured  decimeter 
length  as  a  wave  period  in  seconds.  The  period  thus  determined  was 
compared  to  the  period  measured  by  timing  the  pendulum. 

An  average  of  at  least  10  waves  was  taken  for  the  maximum 
wave  height,  maximum  horizontal  force  and  each  of  the  four 
vertical  maximum  vertical  forces  (one  for  each  supporting  wire). 
The  upward  force  for  each  wire  was  measured  separately  from  its 
downward  force.  The  force  in  the  direction  of  wave  travel  was 
measured  separately  from  the  force  in  the  opposite  direction  of 
wave  travel . 

Once  the  raw  data  was  reduced  to  numerical  values,  it  was 
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reduced  by  use  of  an  IBM  model  360/65  digital  computer.  The  inputs 
to  the  computer  program  were  run  serial  number,  wave  period,  water 
depth,  wave  height,  measured  horizontal  force  in  the  direction  of 
wave  travel,  measured  horizontal  force  in  the  direction  opposite 
from  wave  travel,  measured  vertical  force  down,  vertical  force  up 
and  phase  angle  between  the  wave  crest  and  maximum  forces.  The 
program  averaged  the  wave  period  for  each  group  of  waves  where  the 
water  depth  was  constant  and  the  period  nominally  constant.  From 
the  average  period,  the  water  depth  and  the  wave  height,  the 
program  computed  the  depth  to  length  (d/L)  ratio,  the  length  to 
depth  (L/d)  ratio  and  the  steepness  (H/L). 

The  units  of  the  measured  forces  were  grams.  These  measured 
forces  were  converted  to  dimensionless  form  by  the  program.  The 
dimensionless  force  was  calculated  as  follows: 


F 
F    =  max  (41) 

DIM    A3H (4U 


Y 


d  2 


where 


FDIM  =  Dimensionless  force 

F    =  Maximum  force  measured 
max 

A    =  Significant  linear  dimension  equal  to  height  of  the 

model  squared  divided  by  length  of  model  in 

direction  of  wave  propagation  for  all  models  except 

flat  Dlate 

A  =  H     =  Heiqht  of  model    for  flat  plate 
m  ^  r 
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The  computer  program  which  calculates  wave  parameters  and  dimension- 
less  forces  and  a  sample  of  computer  output  are   presented  in 
Appendix  VI. 

A  dimensionless  plot  was  made  of  dimensionless  force  versus 
steepness  for  each  group  of  waves  with  constant  depth  and  period 
(constant  L/d  ratio).  A  sample  plot  of  dimensionless  force  versus 
steepness  is  shown  in  Fig.  12.  The  horizontal  force  with  the 
direction  of  the  wave  was  usually  slightly  hiaher  than  the  hori- 
zontal force  exerted  against  the  direction  of  the  wave.  As  the 
forces  in  the  two  directions  closely  approximated  each  other  in 
magnitude,  a  single  curve  was  used  to  represent  both  sets  of  points. 

A  cross-plot  was  made  from  the  dimensionless  force  versus 
steepness  plot.  This  cross-plot  summarized  data  by  plotting 
dimensionless  force  versus  L/d  ratio  for  wave  steepnesses  of  0.02  to 
0.07.  These  summary  plots  are   made  for  vertical  and  horizontal 
forces  for  all  models  with  the  exception  of  the  flat  plate.  Only 
horizontal  forces  were  plotted  for  the  flat  plate.  These  curves 
are   presented  in  Appendix  II.  Since  data  are  in  dimensionless 
form,  they  can  be  extrapolated  to  any  size  rectangular  object  or 
to  any  size  of  a  geometrically  similar  half-cylinder  object  for 
any  wave  height,  wave  length  and  water  depth.  An  illustration  of 
the  use  of  the  curves  is  given  in  Appendix  V. 

An  attempt  was  made  to  evaluate  both  inertial  and  drag 
coefficient  by  use  of  the  Morison  Eguation.  This  evaluation  was 
made  by  computer  solution  of  Eq.  28,  Eg.  31  and  Eg.  34.  The 
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inputs  to  this  solution  were  the  same  as  those  cited  above  for 

dimensionless  force  computation.  The  proqram  increased  the 

values  of  C.  and  C  by  iteration.  This  was  done  by  iterating 
dm 

larger  values  of  C, ,  then  holding  them  constant,  and  comparing 
the  total  measured  force  and  the  computed  total  force  for  in- 
creasing values  of  C  .  When  the  measured  force  and  computed  force 
approximate  each  other,  the  comDUted  value  of  C  is  given.  In 
the  case  where  drag  force  is  significant  the  best  values  of  C  ■ 

and  C  can  be  selected  by  picking  the  values  which  remain  more 
m 

or  less  constant  with  changing  wave  parameters.  For  the  flat 
plate  and  two  inch  long  rectangular  model,  the  volume  of  a  cir- 
cumscribing cylinder  was  used  for  the  computations.  It  was  found 
that  drag  forces  were  small  in  the  case  of  the  5/8  inch  thick 
flat  plate.  A  value  of  3.5  for  C,  was  suggested  by  Brater,  McNown 
and  Stair  (1).  Using  this  value,  the  drag  force  represented 
less  than  5  percent  of  the  total  force  for  waves  of  small  steep- 
ness and  10  percent  of  the  total  force  for  waves  of  large  steep- 
ness. Apparently,  the  finite  thickness  of  the  plate  was  enough 
to  make  inertial  force  predominant.  Due  to  the  small  size  of  the 
drag  forces  it  was  impossible  to  evaluate  drag  coefficient.  For 
all  models  except  the  flat  plate,  the  drag  force  was  insignificant, 
The  computer  program  which  applies  the  Mori  son  Equation  to  deter- 
mine drag  and  inertial  coefficients  and  a  sample  of  computer 
output  are  presented  in  Appendix  VIII. 

Since  forces  on  all  the  models  are  largely  inertial,  the 
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values  of  coefficient  of  inertia  were  computed  following  the 

methods  suggested  by  Reid  and  Bretschneider  (15).  A  computer 

program  was  used  which  utilized  the  same  inputs  as  those  cited 

above  for  dimensionless  force  computation.  The  program  determined 

C  by  use  of  Eq.  40.  Values  of  C  were  increased  by  iteration 
m  m 

until  measured  force  and  the  computed  force  aDproximated  each 
other.  This  method  was  used  to  determine  inertial  coefficients 
for  all  models  except  the  flat  plate.  The  computer  program  which 
computes  C  utilizing  Eq.  40  and  a  sample  of  computer  output 
are  presented  in  Appendix  VI.  The  subprogram  titled  "Force-2" 
performs  this  function.  The  program  was  modified  to  sum  a  series 
of  thin  rectangular  solids  to  accomodate  the  half  cylinder  model, 
as  it  was  not  of  constant  cross  section  in  the  direction  of  wave 
travel.  This  program  can  be  applied  to  any  large  submerged  object 
of  irregular  cross  section.  The  program  thus  modified  and  a  samDle 
of  computer  output  are  presented  in  Appendix  VII.  The  subprogram 
"Force-4"  performs  this  function.  The  results  of  calculation  of 
C  by  Eq.  40  are  presented  graphically  in  Appendix  III.  These 
curves  indicate  that  the  coefficients  of  inertia  show  a  small 
variation  with  L/d  and  H/L  ratios.  In  general  the  coefficient 
was  higher  for  higher  values  of  H/L.  There  was  also  some 

variation  of  C  with  L/d  as  indicated  in  the  curves.  The  values 

m 

of  coefficient  C  varied  between  1.20  and  2.20  for  all  tests. 

m 

Methods  suggested  by  Reid  and  Bretschneider  (15)  were  used  to 
compare  the  computed  and  experimental  values.  Several  values  of 
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C  were  selected  for  this  comparison.  The  computation  utilized 
m 

Eq.  40  to  determine  the  calculated  force  for  a  qiven  C  .  Sub- 

m 

program  "Force-1"  performed  this  computation  in  the  computer 
program  for  objects  of  regular  cross  section.  This  program  is 
presented  in  Appendix  VI.  Subprogram  "Force-3"  performed  this 
computation  in  the  computer  program  for  objects  of  irregular  cross 
section.  This  program  is  presented  in  Appendix  VII.  The  results 
of  the  comparison  of  measured  force  and  calculated  force  are 
presented  in  tabular  form  in  Appendix  IV.  The  best  correlation 
of  the  computed  and  experimental  data  for  the  two-inch  long  model 

was  obtained  with  a  value  of  C  of  1.40.  The  best  correlation 

m 

for  all  other  models  was  obtained  with  a  value  of  C  of  1.80. 

m 
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CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

It  was  found  that  force  was  a  regular  function  of  steepness 
(H/L)  and  relative  wave  length  (L/d).  The  force  increased  for 
increasing  values  of  either  parameter.  The  general  shape  was  the 
same  for  all  curves  of  steepness  versus  dimensionless  force. 
The  shape  of  all  relative  wave  length  versus  dimensionless  force 
was  similar.  As  the  data  is  in  dimensionless  form,  it  can  be 
extrapolated  to  prototypes  of  similar  shape. 

The  data  was  analyzed  utilizing  the  Morison  Equation.  It 
was  found  that  drag  force  was  small  for  a  flat  plate  of  finite 
thickness  and  negligible  for  models  with  larger  volumes.  The 
inertial  portion  of  the  Morison  Equation  and  equations  derived 
from  it  were  used  for  analysis  of  forces  and  determination  of 
inertial  coefficient.  The  coefficients  of  inertia  varied  little 
with  changes  in  relative  wave  length  and  steepness.  In  general 
the  coefficient  was  slightly  higher  for  higher  steepnesses.  The 
values  of  inertial  coefficient  varied  between  1.20  and  2.20 
for  all  tests.  Values  of  CM  in  this  range  were  used  to  compute 
force.  This  force  was  compared  to  the  measured  values.  The  best 
correlation  for  the  shortest  (in  the  direction  of  wave  travel) 
model  was  obtained  with  a  value  of  CM  of  1.40.  The  best  cor- 
relation for  all  other  models  was  obtained  with  a  value  of  C..  of 
1.80. 

As  the  flat  plate  tested  yielded  primarily  inertial  forces, 
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future  experiments  should  include  tests  of  a  thinner  plate.  This 
would  facilitate  evaluation  of  drag  coefficient. 

Prototypes  are  exposed  to  the  irregular  gravity  waves  of  the 
sea,  and  not  the  regular  ones  of  the  laboratory.  For  this  reason, 
future  research  should  include  tests  of  models  exposed  to  irregular, 
wind -generated  waves.  These  tests  could  be  related  to  existinq 
work  by  techniques  of  spectral  analysis. 

It  is  also  suggested  that  future  tests  be  conducted  in  a 
three-dimensional  wave  tank  to  eliminate  any  possible  errors 
introduced  by  reflections  induced  by  the  walls  of  a  two-dimensional 
wave  channel . 

The  most  valuable  future  research  would  be  the  instrumentation 
of  prototypes  and  \/ery   large  models.  This  has  already  provided 
valuable  information  in  the  case  of  offshore  platforms.  While 
it  is  a  difficult  and  expensive  proposal,  only  when  prototypes 
are  tested  will  absolute  information  about  forces  on  submerged 
structures  be  known. 
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APPENDIX  I 
NOTATION 


The  following  symbols  are  used  in  this  paper: 

Symbol  Quantity  Dimension 

a         Wave  amplitude  =  H/2  L 

A         Significant  linear  dimension  eoual  to 
heiqht  of  the  model  squared  divided 
by  lenqth  of  model  in  direction  of 
wave  propaqation  for  all  models  except 
flat  olate  L 

A         Height  of  model  for  flat  plate  L 

C  Celerity  L/T 

d         Water  depth  L 

2 

Fr,        Force  due  to  drag  ML/T 

Fpj„      Dimensionless  force 

2 

Fh  Horizontal   force  ML/T 

2 

F,        Force  due  to  inertia  ML/T 

2 
F         Maximum  force  measured  ML/T 

2 
q         Gravity  L/T 

H         Wave  heiqht  L 

H         Heiqht  of  model  L 

m 

K         Pressure  coefficient 


These  symbols  are  consistent  throuoHout  the  thesis  with  the 
exception  of  the  literature  survey.  In  the  literature  survey, 
symbols  are  used  as  they  appear  in  the  cited  literature,  and  are 
redefined  for  each  author. 
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Symbol  Quantity 

k  ?7t/L 

L  Wave  lenqth 

L         Lenath  of  model 
m 

L         Deep  water  wave  lenqth 

P         Gauae  pressure 

AP        Pressure  anomaly 

aP,       Pressure  anomaly  at  leading  edae  of 
object 

AP?       Pressure  anomaly  at  trail inq  edqe  of 
object 

(a  )  Pressure  anomaly  under  the  wave  crest 

(a  ).  Pressure  anomaly  under  the  wave  trouqh 

T  Wave  oeriod 

t  Time 

u  Horizontal  velocity  of  water  oarticle 

u         Horizontal  acceleration  of  water 
particle 

W         Width  of  model 
m 

x         The  distance  from  the  wave  crest 
measured  in  the  direction  of  wave 
travel  L 

x.        Horizontal  distance  from  the  wave  crest 
to  the  edae  of  the  model  direction  of 
wave  travel  L 

z         Hi  stance  from  still  water  level  L 

2  2 

y         Snecific  weiaht  of  water  M/L  T 

n         Distance  of  the  free  surface  -from 

mean  water  level  L 


Dimension 

1/L 

L 

L 

L 

M/LT2 

M/LT2 

M/LT2 

M/LT2 

M/LT2 

M/LT2 

T 

T 

L/T 

2 

in 

L 
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Symbol 


Quantity 


Density 

2tt/T 

Velocity  potential 


ni men si  on 


M/L' 


1/T 

L2/T 


APPFMnix  II 


MMENSIONLFSS   FORCE  SUMMARY  CIJPVFS 
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FIG.I3-DIMENSI0NLESS      HORIZONTAL     FORCE     PLOT     FOR 
FLAT     PLATE. 
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FIG.  I4-DIMENSI0NLESS     HORIZONTAL     FORCE     (COMPUTED     AS  A 
RECTANGULAR  MODEL)   FOR  2  INCH   RECTANGULAR 
MODEL 
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FIG.    I5-DIMENSI0NLESS     VERTICAL    FORCE    (COMPUTED  AS  A 
RECTANGULAR    MODEL)     FOR  2  INCH  RECTANGULAR 
MODEL 
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FIG.  16  -DIMENSIONLESS     HORIZONTAL     FORCE     (COMPUTED     AS 

A    FLAT     PLATE)     FOR     2     INCH     RECTANGULAR     MODEL 
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FIG.  I7-DIMENSI0NLESS     VERTICAL    FORCE    (COMPUTED  AS 
A    FLAT   PLATE)    FOR    2    INCH    RECTANGULAR 
MODEL 
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FIG.  I8-DIMENSI0NLESS     HORIZONTAL     FORCE    FOR  4  INCH 
RECTANGULAR    MODEL 
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FIG.  I9-CHMENSI0NLESS     VERTICAL     FORCE     FOR     4     INCH 
RECTANGULAR      MODEL 
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FIG.  20-DIMENSI0NLESS      HORIZONTAL     FORCE     FOR      8 
INCH      RECTANGULAR      MODEL 
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FIG.  21-  DIMENSIONLESS    VERTICAL  FORCE   FOR 
8    INCH    RECTANGULAR   MODEL 
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FIG.  22- DIMENSIONLESS      HORIZONTAL     FORCE     PLOT      FOR 
HALF      CYLINDER      MODEL     (END     FACING     INCIDENT 
WAVE) 
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FIG.24-DIMENSI0NLESS     HORIZONTAL     FORCE     PLOT     FOR 
HALF     CYLINDER    MODEL    (SIDE     FACING     INCIDENT 
WAVE) 
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FIG.25-DIMENSI0NLESS     VERTICAL      FORCE     PLOT    FOR     HALF 
CYLINDER      MODEL     (SIDE     FACING      INCIDENT    WAVE) 


APPENDIX  III 


TNERTIAL  COEFFICIENT  CIIPVFS 


67 


2.0 


•m 


1.0 


0.0 


*i 

* 

*Vi 

-4,-U 

i 

«•**" 

*•:_. 

u. 

o 



^a, 

n^r&, 

£/ 

SYMBOL 

L/d 
2.66 

— o 

4.52 

— -i 

V — 

5.86 

0.0I   0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09  0.I0   O.I  I   O.I 2 

H/L 
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FIG.  27 -COEFFICIENT     OF     INERTIA     FOR      RECTANGU- 
LAR     MODEL      2 -INCHES      LONG     (WATER 
DEPTH    18    INCHES) 
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FIG.  28- COEFFICIENT     OF     INERTIA     FOR     RECTANGU- 
LAR    MODEL     2  -INCHES      LONG     (WATER 
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FIG.  29- COEFFICIENT     OF     INERTIA     FOR     RECTANGU 
LAR     MODEL     4  -  INCHES      LONG     (WATER 
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FIG.  30- COEFFICIENT    OF    INERTIA     FOR     RECTANGU- 
LAR     MODEL      4 -INCHES      LONG     (WATER 
DEPTH    18   INCHES) 
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FIG.  31- COEFFICIENT     OF     INERTIA     FOR     RECTANGU- 
LAR     MODEL      4- INCHES      LONG     (WATER 
DEPTH  24   INCHES) 
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FIG.    32-   COEFFICIENT    OF    INERTIA    FOR 
RECTANGULAR   MODEL  8-INCHES 
LONG    (WATER    DEPTH    13  INCHES) 
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FIG.  33-  COEFFICIENT  OF  INERTIA  FOR 
RECTANGULAR  MODEL  8-INCHES 
LONG  (WATER    DEPTH    18   INCHES) 
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FIG.  36-COEFFICIENT    OF    INERTIA  FOR  HALF 

CYLINDER   MODEL  (END   FACING  INCIDENT 
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FIG.  39-C0EFICIENT    OF   INERTIA    FOR   HALF 

CYLINDER    MODEL  (SIDE  FACING  INCIDENT 
WAVE,  WATER  DEPTH   18  INCHES) 
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APPENDIX  IV 
COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL  FORCES 


COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES  FOR  RECTANGULAR  MODEL  8  INCHES  LONG 


d  in 

d/L 

L/d 

H/L 

Measured 

Compt 

jted  Horizontal 

feet 

Horizontal 
Force  in 

Foi 

*ce  in  qrams 

grams 

C  =1.5 
m 

C  =1.8 
m 

cm=2-n 
m 

2.000 

0.182 

5.496 

0.007 

94 

89 

107 

119 

2.000 

0.182 

5.496 

0.016 

224 

194 

232 

258 

2.000 

0.182 

5.496 

0.023 

366 

276 

331 

368 

2.000 

0.182 

5.496 

0.032 

528 

387 

465 

517 

2.000 

0.182 

5.496 

0.044 

720 

534 

641 

713 

2.000 

0.182 

5.496 

0.056 

850 

690 

828 

920 

2.000 

0.182 

5.496 

0.065 

1000 

791 

949 

1050 

2.000 

0.354 

2.824 

0.030 

100 

132 

159 

176 

2.000 

0.354 

2.824 

0.045 

210 

201 

242 

286 

2.000 

0.354 

2.824 

0.067 

342 

299 

359 

398 

2.000 

0.354 

2.824 

0.092 

460 

409 

491 

545 

2.000 

0.354 

2.824 

0.092 

524 

407 

488 

542 

2.000 

0.354 

2.824 

0.103 

585 

456 

547 

608 

2.000 

0.656 

1.524 

0.026 

20 

17 

20 

22 

2.000 

0.656 

1.524 

0.051 

42 

33 

39 

43 

2.000 

0.656 

1.524 

0.074 

66 

47 

56 

62 

2.000 

0.656 

1.524 

0.100 

96 

63 

76 

85 

2.000 

0.656 

1.524 

0.112 

106 

71 

85 

95 

2.000 

0.656 

1.524 

0.112 

98 

71 

85 

95 

1.500 

0.577 

1.734 

0.023 

26 

24 

28 

31 

1.500 

0.577 

1.734 

0.045 

56 

45 

54 

61 

1.500 

0.577 

1.734 

0.066 

88 

67 

81 

90 

1.500 

0.577 

1.734 

0.088 

122 

89 

107 

119 

1.500 

0.577 

1.734 

0.111 

146 

113 

136 

151 

1.500 

0.577 

1.734 

0.106 

162 

107 

129 

143 
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COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES  FOR  RECTANGULAR  MODEL  8  INCHES  LONG 
(Continued) 


d  in 

d/L 

L/d 

H/L 

Measured 

Computed  Horizontal 

feet 

Horizontal 
Force  in 

Force  in  grams 

grams 

C  =1.5 

m 

C  =1.8 

m 

Cm=2-0 

m 

1.500 

0.171 

5.851 

0.009 

138 

118 

141 

157 

1.500 

0.171 

5.851 

0.021 

310 

267 

321 

357 

1.500 

0.171 

5.851 

0.031 

516 

400 

480 

533 

1.500 

0.171 

5.851 

0.039 

708 

508 

610 

678 

1.500 

0.171 

5.851 

0.055 

920 

711 

853 

948 

1.500 

0.171 

5.851 

0.063 

1080 

811 

973 

1080 

1.500 

0.171 

5.851 

0.074 

1200 

949 

1140 

1270 

1.083 

0.167 

6.002 

0.012 

160 

153 

183 

204 

1.083 

0.167 

6.002 

0.012 

170 

153 

183 

204 

1.083 

0.167 

6.002 

0.012 

192 

155 

186 

206 

1.083 

0.167 

6.002 

0.025 

408 

322 

386 

429 

1.083 

0.167 

6.002 

0.038 

644 

499 

598 

665 

1.083 

0.167 

6.002 

0.051 

952 

668 

801 

890 

1.083 

0.167 

6.002 

0.069 

1256 

901 

1080 

1200 

1.083 

0.167 

6.002 

0.065 

1152 

849 

1020 

1130 

1 .  083 

0.223 

4.483 

0.022 

234 

213 

256 

284 

1.083 

0.223 

4.483 

0.044 

568 

426 

511 

568 

1.083 

0.223 

4.483 

0.066 

920 

629 

755 

839 

1.083 

0.223 

4.483 

0.091 

1200 

876 

1050 

1170 

1.083 

0.223 

4.483 

0.068 

900 

655 

786 

873 

1.083 

0.424 

2.360 

0.038 

140 

99 

119 

132 

1.083 

0.424 

2.360 

0.075 

288 

198 

238 

264 

1.083 

0.424 

2.360 

0.099 

384 

261 

313 

348 

1.083 

0.424 

2.360 

0.099 

424 

260 

312 

346 

1.083 

0.424 

2.360 

0.091 

364 

239 

287 

319 

1.083 

0.424 

2.360 

0.088 

676 

231 

277 

308 

COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES    FOR   RECTANGULAR  MODEL   4   INCHES  LONG 
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d  in 

d/L 

'   L/d 

H/L 

Measured 

Computed  Horizi 

Dntal 

feet 

Horizontal 

Force  in  grams 

Force  in 
qrams 

C  =1.6 

C  =1.8 

C  =2.0 

j 

m 

m 

m 

2.000 

0.173 

5.770 

0.031 

320 

323 

363 

404 

2.000 

0.173 

5.770 

0.038 

431 

389 

437 

486 

2.000 

0.173 

5.770 

0.045 

520 

468 

526 

584 

2.000 

0.173 

5.770 

0.045 

230 

236 

266 

295 

2.000 

0.173 

5.770 

0.014 

160 

149 

168 

186 

2.000 

0.345 

2.898 

0.037 

180 

141 

158 

176 

2.000 

0.345 

2.898 

0.055 

280 

210 

237 

263 

2.000 

0.345 

2.898 

0.074 

374 

283 

318 

354 

2.000 

0.345 

2.898 

0.075 

373 

287 

323 

359 

2.000 

0.345 

2.898 

0.068 

315 

259 

292 

324 

1.500 

0.200 

4.992 

0.013 

110 

120 

135 

150 

1.500 

0.200 

4.992 

0.032 

245 

287 

323 

359 

1.500 

0.200 

4.992 

0.042 

375 

372 

419 

465 

1.500 

0.200 

4.992 

0.059 

554 

527 

593 

658 

1.500 

0.200 

4.992 

0.067 

716 

601 

676 

751 

1.500 

0.200 

4.992 

0.082 

800 

733 

824 

916 

1.500 

0.479 

2.086 

0.033 

64 

54 

60 

67 

1.500 

0.479 

2.086 

0.066 

132 

107 

121 

134 

1.500 

0.479 

2.086 

0.096 

200 

155 

174 

194 

1.500 

0.479 

2.086 

0.100 

210 

162 

182 

202 

1.500 

0.479 

2.086 

0.094 

200 

152 

171 

190 

1.500 

0.479 

2.086 

0.087 

179 

141 

159 

176 

1.083 

0.378 

2.642 

0.032 

112 

97 

109 

121 

1.083 

0.378 

2.642 

0.067 

230 

201 

226 

252 

1.083 

0.378 

2.642 

0.070 

230 

210 

236 

262 

1.083 

0.378 

2.642 

0.098 

320 

294 

330 

367 

1.083 

0.378 

2.642 

0.098 

330 

294 

330 

367 

1.083 

0.378 

2.642 

0.091 

300 

273 

307 

341 

1.083 

0.378 

2.642 

0.089 

310 

268 

302 

336 

1.083 

0.378 

2.642 

0.085 

304 

255 

287 

318 

85 


COMPARISON  OF     MEASURED  AND  COMPUTED  HORIZONTAL 

FORCES   FOR   RECTANGULAR  MODEL   4  INCHES   LONG 

(Continued) 


d  in 

d/L 

L/d 

H/L 

Measured 

Computed  Horizontal 

feet 

Horizontal 

Force   in  grams 

Force  in 
grams 

C  =1.6 

C  =1.8 

C  =2.0 

m 

m 

m 

1.083 

0.220 

4.543 

0.020 

140 

162 

182 

202 

1.083 

0.220 

4.543 

0.039 

300 

308 

346 

385 

1.083 

0.220 

4.543 

0.057 

460 

454 

510 

567 

1.083 

0.220 

4.543 

0.071 

600 

567 

638 

709 

1.083 

0.220 

4.543 

0.073 

625 

580 

652 

725 

1.083 

0.220 

4.543 

0.086 

750 

685 

771 

856 

1.083 

0.169 

5.910 

0.011 

124 

118 

133 

148 

1.083 

0.169 

5.910 

0.024 

265 

249 

280 

312 

1.083 

0.169 

5.910 

0.039 

457 

412 

463 

514 

1.083 

0.169 

5.910 

0.055 

694 

574 

646 

717 

1  .083 

0.169 

5.910 

0.069 

828 

720 

810 

900 

1.083 

0.169 

6.910 

0.066 

806 

697 

784 

871 

86 


COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES  FOR  RECTANGULAR  MODEL  2  INCHES  LONG 


d  in 

d/L 

L/d 

H/L 

Measured 

Comp 

>uted  Horizontal 

feet 

Horizontal 
Force  in 

Force  in  grams 

grams 

C  -1.2 
m 

C  =1 .4 

m 

Cm=1'6 

m 

2.000 

0.180 

5.547 

0.006 

56 

49 

57 

65 

2.000 

0.180 

5.547 

0.014 

116 

108 

126 

144 

2.000 

0.180 

5.547 

0.015 

118 

112 

131 

150 

2.000 

0.180 

5.547 

0.025 

196 

184 

214 

245 

2.000 

0.180 

5.547 

0.035 

296 

262 

306 

350 

2.000 

0.180 

5.547 

0.043 

392 

320 

374 

427 

2.000 

0.180 

5.547 

0.060 

506 

447 

522 

596 

2.000 

0.180 

5.547 

0.061 

624 

460 

536 

613 

2.000 

0.344 

2.909 

0.021 

58 

61 

71 

81 

2.000 

0.344 

2.909 

0.043 

121 

125 

146 

167 

2.000 

0.344 

2.909 

0.068 

183 

197 

230 

263 

2.000 

0.344 

2.909 

0.089 

244 

259 

302 

345 

2.000 

0.344 

2.909 

0.094 

296 

273 

318 

364 

2.000 

0.344 

2.909 

0.094 

266 

272 

317 

362 

2.000 

0.646 

1.547 

0.024 

5 

10 

12 

14 

2.000 

0.646 

1.547 

0.049 

23 

21 

24 

28 

2.000 

0.646 

1.547 

0.071 

35 

30 

35 

40 

2.000 

0.646 

1.547 

0.094 

48 

40 

47 

54 

2.000 

0.646 

1.547 

0.110 

52 

47 

54 

62 

2.000 

0.646 

1.547 

0.109 

58 

46 

54 

62 

1.500 

0.467 

2.140 

0.032 

43 

42 

50 

57 

1.500 

0.467 

2.140 

0.062 

94 

82 

95 

109 

1.500 

0.467 

2.140 

0.083 

132 

109 

127 

145 

1.500 

0.467 

2.140 

0.099 

142 

130 

152 

173 

1.500 

0.467 

2.140 

0.088 

no 

116 

135 

154 

1.500 

0.201 

4.975 

0.013 

78 

84 

98 

113 

1.500 

0.201 

4.975 

0.027 

182 

178 

207 

237 

1.500 

0.201 

4.975 

0.040 

290 

270 

315 

360 

1.500 

0.201 

4.975 

0.054 

462 

364 

425 

486 

1.500 

0.201 

4.975 

0.061 

568 

411 

480 

548 

1.500 

0.201 

4.975 

0.078 

640 

524 

612 

699 

1.500 

0.201 

4.975 

0.069 

524 

460 

536 

613 

87 


COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES  FOR  RECTANGULAR  MODEL  2  INCHES  LONG 
(Continued) 


d  in 

d/L 

L/d 

H/L 

Measured 

Computed  Horizon 

tal 

feet 

Horizontal 
Force  in 

Force 

»  in  qrams 

C  =1.2 

C  =1-4 

Cm=l .6 

grams 

m 

m 

m 

1.083 

0.171 

5.862 

0.011 

82 

86 

101 

115 

1.083 

0.171 

5.862 

0.022 

184 

174 

203 

231 

1.083 

0.171 

5.862 

0.037 

364 

293 

342 

391 

0.083 

0.171 

5.862 

0.040 

392 

309 

361 

412 

1.083 

0.171 

5.862 

0.057 

564 

446 

520 

594 

1.083 

0.171 

5.862 

0.070 

692 

549 

641 

732 

1.083 

0.171 

5.862 

0.068 

684 

529 

618 

706 

1.083 

0.171 

5.862 

0.064 

644 

500 

583 

667 

1.083 

0.221 

4.523 

0.019 

118 

112 

130 

149 

1.083 

0.221 

4.523 

0.043 

258 

257 

300 

343 

1.083 

0.221 

4.523 

0.067 

422 

397 

463 

529 

1.083 

0.221 

4.523 

0.093 

572 

556 

648 

741 

1.083 

0.221 

4.523 

0.071 

440 

422 

493 

563 

1.083 

0.377 

2.655 

0.031 

78 

70 

81 

93 

1.083 

0.377 

2.655 

0.059 

152 

135 

157 

179 

1.083 

0.377 

2.655 

0.093 

220 

211 

246 

281 

1.083 

0.377 

2.655  - 

0.094 

222 

214 

249 

285 

1.083 

0.377 

2.655 

0.077 

174 

175 

204 

233 

88 


COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES    FOR  HALF  CYLINDER  MODEL    (END   FACING   INCIDENT  WAVE) 


d  in 

d/L 

L/d 

H/L 

Measured 

Comput 

sd  Horizontal 

feet 

Horizontal 
Force  in 

Force  in  arams 

grams 

C  =1.6 

m 

C  =1.8 

m 

Cm=2'0 

m 

2.000 

0.180 

5.552 

0.007 

82 

83 

93 

104 

2.000 

0.180 

5.552 

0.016 

206 

179 

202 

224 

2.000 

0.180 

5.552 

0.024 

344 

277 

311 

346 

2.000 

0.180 

5.552 

0.034 

510 

386 

434 

482 

2.000 

0.180 

5.552 

0.047 

720 

545 

613 

682 

2.000 

0.180 

5.552 

0.055 

800 

638 

718 

798 

2.000 

0.180 

5.552 

0.069 

988 

790 

889 

987 

2.000 

0.351 

2.850 

0.020 

88 

84 

94 

105 

2.000 

0.351 

2.850 

0.042 

184 

179 

202 

224 

2.000 

0.351 

2.850 

0.064 

286 

270 

304 

337 

2.000 

0.351 

2.850 

0.084 

398 

355 

399 

443 

2.000 

0.351 

2.850 

0.097 

516 

411 

462 

514 

2.000 

0.351 

2.850 

0.091 

448 

383 

431 

478 

2.000 

0.653 

1  .532 

0.021 

20 

13 

14 

16 

2.000 

0.653 

1.532 

0.047 

36 

29 

32 

36 

2.000 

0.653 

1.532 

0.070 

56 

43 

48 

54 

2.000 

0.653 

1.532 

0.093 

76 

57 

64 

71 

2.000 

0.653 

1.532 

0.180 

86 

66 

74 

82 

2.000 

0.653 

1.532 

0.100 

86 

61 

68 

76 

1.500 

0.474 

2.108 

0.033 

72 

62 

69 

77 

1.500 

0.474 

2.108 

0.067 

141 

125 

141 

156 

1.500 

0.474 

2.108 

0.092 

212 

173 

194 

216 

1.500 

0.474 

2.108 

0.098 

220 

183 

206 

229 

1.500 

0.474 

2.108 

0.087 

162 

162 

183 

203 

1.500 

0.200 

5.006 

0.011 

127 

113 

127 

141 

1.500 

0.200 

5.006 

0.025 

288 

254 

286 

318 

1.500 

0.200 

5.006 

0.040 

480 

416 

468 

520 

1.500 

0.200 

5.006 

0.055 

724 

566 

637 

708 

1.500 

0.200 

5.006 

0.070 

984 

717 

807 

897 

1.500 

0.200 

5.006 

0.080 

1072 

827 

931 

1030 

1.500 

0.200 

5.006 

0.068 

864 

698 

785 

873 

89 


COMPARISON  OF  MEASURED  AND   COMPUTED  HORIZONTAL 
FORCES   FOR  HALF  CYLINDER  MODEL    (END   FACING   INCIDENT   WAVE) 

(Conti  nued) 


d  in 

d/L 

L/d 

H/L 

Measured 

Computed  Horizontal 

feet 

Horizontal 
Force  in 

Force  in  grams 

grams 

Cm=1'6 

Cm=1'8 

C  =2.0 

m 

1.083 

0.169 

5.926 

0.012 

128 

150 

169 

188 

1  .083 

0.169 

5.926 

0.024 

320 

285 

321 

357 

1.083 

0.169 

5.926 

0.038 

560 

464 

522 

580 

1.083 

0.169 

5.926 

0.059 

904 

714 

803 

892 

1.083 

0.169 

5.926 

0.069 

1152 

826 

930 

1030 

1.083 

0.169 

5.926 

0.062 

1096 

751 

845 

939 

1.083 

0.222 

4.512 

0.018 

172 

162 

183 

203 

1.083 

0.222 

4.512 

0.041 

406 

369 

415 

461 

1.083 

0.222 

4.512 

0.064 

684 

581 

653 

726 

1.083 

0.222 

4.512 

0.088 

952 

797 

896 

996 

1.083 

0.222 

4.512 

0.067 

720 

601 

676 

751 

1.083 

0.382 

2.619 

0.032 

104 

104 

117 

130 

1.083 

0.382 

2.619 

0.065 

220 

212 

239 

266 

1.083 

0.382 

2.619 

0.096 

336 

314 

353 

393 

1.083 

0.382 

2.619 

0.091 

300 

297 

334 

371 

1.083 

0.382 

2.619 

0.086 

288 

282 

317 

352 

90 


COMPARISON  OF  MEASURED  AND  COMPUTED  HORIZONTAL 
FORCES  FOR  HALF  CYLINDER  MODEL  (SIDE  FACING  INCIDENT  WAVE) 


d  in 

d/L 

L/d 

H/L 

Measured 

Compi 

ited  Horizontal 

feet 

Horizontal 

Foi 

xe   in  grams 

Force   in 

grams 

C  =1.6 
m 

Cm=1-8 

C  =2.0 

m 

2.000 

0.179 

5.578 

0.008 

80 

89 

100 

111 

2.000 

0.179 

5.578 

0.016 

204 

179 

202 

224 

2.000 

0.179 

5.578 

0.025 

340 

284 

320 

356 

2.000 

0.179 

5.578 

0.034 

460 

390 

438 

487 

2.000 

0.179 

5.578 

0.045 

600 

515 

580 

644 

2.000 

0.179 

5.578 

0.056 

728 

643 

723 

804 

2.000 

0.179 

5.578 

0.066 

892 

754 

849 

943 

2.000 

0.344 

2.907 

0.021 

78 

91 

102 

113 

2.000 

0.344 

2.907 

0.043 

170 

189 

212 

236 

2.000 

0.344 

2.907 

0.065 

272 

286 

322 

358 

2.000 

0.344 

2.907 

0.087 

380 

382 

430 

477 

2.000 

0.344 

2.907 

0.086 

400 

377 

425 

472 

2.000 

0.344 

2. 907 

0.096 

484 

423 

475 

528 

2.000 

0.344 

2.907 

0.095 

448 

417 

469 

521 

2.000 

0.653 

1.532 

0.021 

18 

13 

14 

16 

2.000 

0.653 

1.532 

0.021 

34 

27 

30 

34 

2.000 

0.653 

1.532 

0.021 

49 

38 

43 

47 

2.000 

0.653 

1.532 

0.021 

69 

50 

56 

62 

2.000 

0.653 

1.532 

0.021 

74 

58 

65 

72 

2.000 

0.653 

1.532 

0.021 

68 

55 

62 

69 

1.500 

0.473 

2.113 

0.030 

61 

57 

64 

71 

1.500 

0.473 

2.113 

0.030 

123 

118 

133 

148 

1.500 

0.473 

2.113 

0.030 

224 

173 

195 

217 

1.500 

0.473 

2.113 

0.030 

238 

169 

190 

212 

1.500 

0.473 

2.113 

0.030 

234 

173 

195 

217 

1.500 

0.473 

2.113 

0.030 

170 

168 

189 

210 

1.500 

0.192 

5.206 

0.010 

121 

110 

123 

137 

1.500 

0.192 

5.206 

0.024 

318 

255 

287 

319 

1.500 

0.192 

5.206 

0.037 

492 

400 

450 

500 

1.500 

0.192 

5.206 

0.050 

720 

537 

604 

671 

1.500 

0.192 

5.206 

0.067 

900 

721 

811 

901 

1.500 

0.192 

5.206 

0.078 

1000 

830 

934 

1004 

91 


COMPARISON  OF  MEASURED  AND  COMPUTED   HORIZONTAL 
FORCES   FOR  HALF  CYLINDER  MODEL   (SIDE   FACING  INCIDENT  WAVE) 

(Continued) 


d  in 

d/L 

L/d 

H/L 

Measured 

Compi 

ited  Horizontal 

feet 

Horizontal 
Force  in 

Force  in  grams 

grams 

C  =1.6 
m 

C  =1.8 

m 

C   =2.0 

m 

1.083 

0.164 

6.088 

0.010 

151 

123 

138 

153 

1.083 

0.164 

6.088 

0.021 

370 

259 

291 

323 

1.083 

0.164 

6.088 

0.037 

604 

450 

506 

563 

1.083 

0.164 

6.088 

0.051 

860 

623 

701 

779 

1.083 

0.164 

6.088 

0.066 

964 

811 

912 

1010 

1.083 

0.164 

6.088 

0.061 

1008 

748 

841 

935 

1.083 

0.164 

6.088 

0.057 

944 

703 

791 

879 

1.083 

0.207 

4.836 

0.017 

174 

164 

184 

205 

1.083 

0.207 

4.836 

0.040 

428 

393 

442 

491 

1.083 

0.207 

4.836 

0.061 

672 

600 

675 

750 

1.083 

0.207 

4.836 

0.080 

916 

779 

876 

973 

1.083 

0.207 

4.836 

0.069 

700 

671 

754 

838 

1.083 

0.367 

2.725 

0.027 

100 

97 

110 

122 

1.083 

0.367 

2.725 

0.063 

260 

226 

255 

283 

1.083 

0.367 

2.725 

0.093 

406 

333 

375 

417 

1.083 

0.367 

2.725 

0.085 

402 

307 

345 

383 

92 


APPENDIX  V 

SAMPLE  COMPUTATION  OF  HORIZONTAL  FORCE 

The  following  is  a  sample  computation  for  a  submerqed  tank 
utilizing  the  di mens ionl ess  force  plots  presented  in  Appendix  II. 
A  check  on  the  answer  obtained  by  the  use  of  these  curves  is  made 
by  computing  the  force  utilizing  the  method  set  forth  by  Reid 
and  Bretschneider  (15). 
Gi  ven : 

Wave:    Period  =  T  =  10  seconds 
Height  =  H  =  10  feet 
Water  depth  =  d  =  150  feet 

Model :         Height  =  H     =  40  feet 

^     m 

Width  =  W  =  60  feet 

m 

Length  =  L  =  70  feet 

3  m 

LQ  =  5.12  T2 (23) 

=  512  feet 

d/LQ  =   150/512  =  0.293 

d/L  =  0.3058  [from  Appendix  I  of  Wiegel  (16)] 

L   =  491 

L/d  =  3.27 

H/L  =  10/491  =  0.02 

^DIM  =  14'3  [from  Appendix  II,  Fig.  20  ] 
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max 


DIM 

Y  A3  H 

d     2 

4.3  = 

max 

Y  A3  H 
d     2 

A.  V     .   (40H40)   .  228 

m 

A3  =   11,850 

(FDIM)(A3)(Y)(H) 
Fmax  JJW) 

=    (14.3)(11850)(62.4)(10) 

mum 

=  352,000  LB 
Check: 

An  inertial    coefficient  of  1.8  will    be  used   in  Eq.   40, 

ttL 


m> 


(Fu)  =  C     H     W     K  Y  H  sin    (-j-^) 

v   Ir  max         m     m     m       '  v   L   ' 


(40) 


K  =  0.30     (from  Fig.  41  or  from  Eq.  37] 

(FH)max  -  (1.8)(40)(60)(0.30)(62.4)(10)(sin(^-)) 

=  (1.8) (40) (60) (0.30) (62. 4) (10) (sin (0.448)) 
=  347,000  LBS 
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APPENDIX  VI 
COMPUTER  PROGRAM  FOR  WAVE  PARAMETERS  AMP  IMERTIAL  WAVE  FORCES 
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APPENDIX  VII 
COMPUTER  PROGRAM  EOR  OBJECTS  OF  IRREGULAR  CROSS  SECTION 
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APPENDIX  VIII 
COMPUTER  PROGRAM  FOR  FORCE  COMPUTATION  RY  MORISON  EQUATION 
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VITA 

George  Edward  Shank  was  born  in  Washington,  D.  C.  on  May  12, 
1937  of  Ralph  and  Evelyn  Shank.  He  received  a  NROTC  Scholar- 
ship and  entered  Duke  University  in  September,  1955.  He  completed 
his  undergraduate  work  in  May  1960  and  received  a  Bachelor  of 
Science  in  Civil  Engineering.  Upon  graduation  he  was  commissioned 
in  the  Civil  Engineer  Corps  of  the  U.  S.  Navy  as  Ensign  and  has 
remained  on  active  duty  since  that  time.  Positions  which  he  has 
held  with  the  Navy  have  been  Project  Manager  for  Aeroballi sties  at 
the  Naval  Ordnance  Laboratory;  Assistant  Resident  Officer  in  Charge 
of  Construction,  Naval  Air  Station,  Memphis,  Tennessee;  Engineer- 
ing Coordination  Officer,  Southeast  Division,  Naval  Facilities 
Engineering  Command;  Public  Works  Officer,  Atlantic  Undersea  Test 
and  Evaluation  Center;  and  Real  Estate  Officer,  Military  Assistant 
Command,  Vietnam. 

In  June  of  1968  he  entered  Texas  A&M  University  to  continue 
his  studies  in  Civil  Engineering.  He  completed  his  studies  and 
was  awarded  a  Master  of  Science  Degree  in  May  1970. 

His  permanent  address  is  4707  Calvert  Road,  College  Park, 
Maryland,  20740. 
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